Method for multi-center effect compensation based on PET/CT intelligent diagnosis system

ABSTRACT

Disclosed is a method for multi-center effect compensation based on a PET/CT intelligent diagnosis system. The method includes the following steps: estimating multi-center effect parameters of a test center B relative to a training center A by implementing a nonparametric mathematical method for data of the training center A and the test center B based on a location-scale model about additive and multiplicative multi-center effect parameters, and using the parameters to compensate the data of the test center B to eliminate a multi-center effect between the test center B and the training center A. According to the present disclosure, the multi-center effect between the training center A and the test center B can be compensated, so that the compensated data of the test center B can be used in the model trained by the training center A, and the generalization ability of the model is indirectly improved.

CROSS-REFERENCE TO RELATED APPLICATIONS

The present application is a National Phase of International Application No. PCT/CN2021/073464, filed on Jan. 23, 2021, which claims priority to Chinese Application No. 202010109082.9, filed on Feb. 21, 2020, the contents of both of which are incorporated herein by reference in their entireties.

TECHNICAL FIELD

The present disclosure relates to the field of medical imaging and deep learning, in particular to a method for multi-center effect compensation based on a PET/CT intelligent diagnosis system.

BACKGROUND

Positron emission tomography (PET) is a functional imaging device at the molecular level. Before scanning, it is necessary to inject a radioactive tracer into a patient. The tracer decays in the patient's body and then annihilates, producing a pair of 511 keV gamma photons with opposite emission directions of about 180°. A detector will collect the position and time information of these gamma photons reaching the crystal. By using an image reconstruction algorithm to reconstruct and post-process the collected information, the metabolism and uptake of the radioactive tracer in the patient can be obtained. According to the image results of PET/CT and various clinical indicators, doctors comprehensively analyze the patient's condition, so as to determine a treatment plan.

A convolutional neural network (CNN) is one of the commonly used methods to build medical artificial intelligence models in recent years. It extracts high-order feature information of images through multi-layer convolution processing, and combines pool processing to reduce the dimension of features. The extracted high-order features are input to a subsequent specific network for specific tasks, such as classification, segmentation, registration, detection, noise reduction and so on. The advantage of this method is that it can automatically learn high-order features of significant significance for specific tasks through a large number of samples, but it has certain requirements for the amount of data used for training.

As for an existing intelligent diagnosis model, whether the model is based on single-center data or multi-center data training, the diagnostic efficiency of other center data that did not participate in the training will decrease for the center data that participate in the training. This is because different centers usually use different scanning protocols, which leads to different levels of image data in different centers. For example, on PET/CT, because different centers use different scanning times per bed, image reconstruction algorithms, iteration times/subset numbers, post-reconstruction filters, voxel sizes and other parameters, the PET/CT image data of these centers will have certain differences, that is, the above multi-center effect. This problem can be solved by fine tuning by adding the data that do not participate in the training center into the model, but this solution has certain requirements for the amount of data. However, the number of cases in grass-roots hospitals often falls far short of the sample size requirement for fine-tuning, so this method is not realistic for these grass-roots hospitals. Therefore, it is an important demand to solve the present disclosure problem of the multi-center intelligent diagnosis model in grass-roots hospitals.

SUMMARY

In view of the shortcomings of the prior art, the aim of the present disclosure is to provide a method for multi-center effect compensation based on a PET/CT intelligent diagnosis system. Using the location-scale model (L/S model) and the empirical Bayesian method, the feature maps of deep neural networks with different centers are standardized to eliminate the multi-center effect encountered in the application of an intelligent diagnosis model. Only the features of all data participating in model training need to be saved, and the features of data from a test center B are used for calculation, so that the multi-center effect parameters between the features of data from a training center A and the test center B can be estimated, and finally, the features of data from the test center B are compensated to the level of the features of data from the center A, and the compensated features of data from the test center B are applied to the intelligent diagnosis system of the training center A.

The purpose of the present disclosure is realized by the following technical solution: a method for multi-center effect compensation based on a PET/CT intelligent diagnosis system, specifically including the following steps:

Step 1, importing PET/CT data of a training center A into a deep convolution neural network model of the training center A, obtaining a feature map of all the PET/CT data in the training center A, calculating an average value {circumflex over (α)}_(g), a covariant regression coefficient {circumflex over (β)}_(g) and a variance {circumflex over (σ)}_(g) of all features Y_(Ajg) in the feature map of all the PET/CT data in the training center A, and carrying out Z-Score standardization thereon to obtain a standardized feature Z_(Ajg):

$\begin{matrix} {Z_{Ajg} = \frac{Y_{Ajg} - {\overset{\hat{}}{\alpha}}_{g} - {X{\overset{\hat{}}{\beta}}_{g}}}{{\overset{\hat{}}{\sigma}}_{g}}} & (1) \end{matrix}$

where Z_(Ajg) represents a feature g extracted from a case j in the training center A, and X is a design matrix of a relevant covariates;

Step 2, importing PET/CT data in a test center B into the deep convolution neural network model of the training center A, obtaining a feature map of all the PET/CT data in the test center B, and carrying out Z-Score standardization on the features Y_(Bjg) in the feature map by using the average value and the variance of the training center A obtained in step 1, and obtaining a standardized feature Z_(Bjg):

$\begin{matrix} {{Z_{Bjg} = \frac{Y_{Bjg} - {\overset{\hat{}}{\alpha}}_{g} - {X{\overset{\hat{}}{\beta}}_{g}}}{{\overset{\hat{}}{\sigma}}_{g}}};} & (2) \end{matrix}$

Step 3, estimating multi-center effect parameters γ*_(ig) and δ*_(ig) according to the standardized features of the training center A and the test center B, and carrying out multi-center effect compensation for the feature of the test center B, so as to obtain a feature Y*_(ijg) of the test center B after the multi-center effect compensation, which comprise that following sub-steps:

(3.1) the standardized feature Z_(Ajg) obtained in step 1 and the standardized feature Z_(Bjg) obtained in step 2 conforming to a normal distribution, and calculating an average value γ_(ig) and a standard deviation {circumflex over (δ)}_(ig) ² of the standardized features; obtaining an expected value E[γ_(ig)] about a multi-center effect parameter γ_(ig) by integration from a posterior distribution π(Z_(ig), γ_(ig), δ_(ig) ²) of standardized data Z_(ig); transforming the integral by a likelihood function according to the Bayesian theorem;

$\begin{matrix} {{\overset{\hat{}}{\gamma}}_{ig} = {\frac{1}{n_{i}}{\sum_{j}Z_{ijg}}}} & (3) \end{matrix}$ $\begin{matrix} {{\overset{\hat{}}{\delta}}_{ig}^{2} = {\frac{1}{n_{i - 1}}{\sum_{j}\left( {Z_{ijg} - {\overset{\hat{}}{\gamma}}_{ig}} \right)^{2}}}} & (4) \end{matrix}$ $\begin{matrix} {{E\left\lbrack \gamma_{ig} \right\rbrack} = {\int{\gamma_{ig}{\pi\left( {Z_{ig},\gamma_{ig},\delta_{ig}^{2}} \right)}{d\left( {\gamma_{ig},\delta_{ig}^{2}} \right)}}}} & (5) \end{matrix}$ $\begin{matrix} {{E\left\lbrack \gamma_{ig} \right\rbrack} = {\frac{1}{C\left( Z_{ig} \right)}{\int{\gamma_{ig}{L\left( {Z_{ig}{❘{\gamma_{ig},\delta_{ig}^{2}}}} \right)}{\pi\left( {\gamma_{ig},\delta_{ig}^{2}} \right)}{d\left( {\gamma_{ig},\delta_{ig}^{2}} \right)}}}}} & (6) \end{matrix}$

where, n_(i) represents a sample size of a center i; Z₁₉ represents the standardized features, in which i represents a center, j represents a case, and g represents a feature, C(Z_(ig))=∫L(Z_(ig)|γ_(ig), δ_(ig) ²)π(γ_(ig), δ_(ig) ²)d(γ_(ig), δ_(ig) ²) L(Z_(ig)|γ_(ig), δ_(ig) ²)=Π_(j)φ(Z_(ijg), γ_(ig), δ_(ig) ²), φ(Z_(ijg), γ_(ig), δ_(ig) ²) is a probability density function for evaluating a certain feature Z_(ijg) about a random variable N(γ_(ig), δ_(ig) ²) that conforms to the normal distribution;

(3.2) calculating an integral about the multi-center effect parameters by a Monte Carlo integral method, defining w_(ig″)=L(Z_(ig)|γ_(ig″), δ_(ig″)), so that an additive multi-center effect parameter γ*_(ig) can be estimated; using a same process to estimate a multiplicative multi-center effect δ*_(ig);

$\begin{matrix} {\gamma_{ig}^{*} = {{\hat{E}\left( \gamma_{ig} \right)} = \frac{\sum_{g^{''}}{w_{{ig}^{''}}{\overset{\hat{}}{\gamma}}_{{ig}^{''}}^{2}}}{\sum_{g^{''}}w_{{ig}^{''}}}}} & (7) \end{matrix}$ $\begin{matrix} {\delta_{ig}^{*} = {{\hat{E}\left( \delta_{ig} \right)} = \frac{\sum_{g^{''}}{w_{{ig}^{''}}{\overset{\hat{}}{\delta}}_{{ig}^{''}}^{2}}}{\sum_{g^{''}}w_{{ig}^{''}}}}} & (8) \end{matrix}$

(3.3) substituting the multi-center effect parameters γ*_(Bg) and δ*_(Bg) estimated in step (3.2) into the location-scale model to obtain the feature Y*_(Bjg) of the test center B after the multi-center effect compensation;

$\begin{matrix} {{Y_{{Bjg}^{=}}^{*}\frac{{\overset{\hat{}}{\sigma}}_{g}}{{\hat{\delta}}_{Bg}^{*}}\left( {Z_{Bjg} - {\overset{\hat{}}{\gamma}}_{Bg}^{*}} \right)} + {\overset{\hat{}}{\alpha}}_{g} + {X{{\overset{\hat{}}{\beta}}_{g}.}}} & (9) \end{matrix}$

The method has the beneficial effects that a location-scale model is used and the empirical Bayesian method is utilized to estimate the multi-center effect parameters, so as to eliminate the multi-center effect among different central data. For a test center B, the amount of data is not enough to fine-tune a deep convolution neural network model, which leads to the result that the model cannot achieve the desired effect in the center B. According to the present disclosure, the multi-center effect between the training center A and the test center B can be reduced, so that the standardized data of the test center B can be used in the model trained by the training center A, and the generalization ability of the model is indirectly improved. This method can indirectly improve the generalization ability of the intelligent diagnosis model, so that the intelligent diagnosis model trained in large-scale top-three hospitals can also be applied in grass-roots hospitals, which is very helpful for improving the diagnosis level of grass-roots hospitals and reducing the workload of doctors.

BRIEF DESCRIPTION OF DRAWINGS

FIG. 1 is a flow chart of the multi-center compensation method of the present disclosure;

FIG. 2 is a flow chart showing the estimation of multi-center effect parameters.

DESCRIPTION OF EMBODIMENTS

The present disclosure will be described in detail below with reference to the attached drawings.

FIG. 1 is a flow chart of a method for multi-center effect compensation based on a PET/CT intelligent diagnosis system of the present disclosure, which specifically includes the following steps:

Step 1, importing PET/CT data of a training center A into a deep convolution neural network model of the training center A, obtaining a feature map of all the PET/CT data in the training center A, calculating an average value {circumflex over (α)}_(g), a covariant regression coefficient {circumflex over (β)}_(g) and a variance {circumflex over (σ)}_(g) of all features Y_(Ajg) in the feature map of all the PET/CT data in the training center A, and carrying out Z-Score standardization thereon to obtain a standardized feature Z_(Ajg):

$\begin{matrix} {Z_{Ajg} = \frac{Y_{Ajg} - {\overset{\hat{}}{\alpha}}_{9} - {X{\overset{\hat{}}{\beta}}_{g}}}{{\overset{\hat{}}{\sigma}}_{g}}} & (1) \end{matrix}$

where Z_(Ajg) represents a feature g extracted from a case j in the training center A, and X is a design matrix of a relevant covariates;

Step 2, importing PET/CT data in a test center B into the deep convolution neural network model of the training center A, obtaining a feature map of all the PET/CT data in the test center B, and carrying out Z-Score standardization on the features Y_(Bjg) in the feature map by using the average value and the variance of the training center A obtained in step 1, and obtaining a standardized feature Z_(Bjg):

$\begin{matrix} {{Z_{Bjg} = \frac{Y_{Bjg} - {\overset{\hat{}}{\alpha}}_{g} - {X{\overset{\hat{}}{\beta}}_{g}}}{{\overset{\hat{}}{\sigma}}_{g}}};} & (2) \end{matrix}$

Because the deep convolutional neural network model is trained based on the data of the training center A, the features of the test center B can be compensated to the level of the training center A through the above two steps.

Step 3, estimating multi-center effect parameters γ*_(ig) and δ*_(ig) according to the standardized features of the training center A and the test center B, and carrying out multi-center effect compensation for the feature of the test center B, so as to obtain a feature Y*_(ijg) of the test center B after the multi-center effect compensation, which comprise that following sub-steps:

(3.1) the standardized feature Z_(Ajg) obtained in step 1 and the standardized feature Z_(Bjg) obtained in step 2 conforming to a normal distribution, and calculating an average value {circumflex over (γ)}_(ig) and a standard deviation δ_(ig) ² of the standardized features to be used as an observation value of the multi-center effect parameter; obtaining an expected value E[γ_(ig)] about a multi-center effect parameter γ_(ig) by integration from a posterior distribution π(Z_(ig), γ_(ig), δ_(ig) ²) of standardized data Z_(ig) (formula 5); substituting the integral of formula 5 into the likelihood function to transform into formula 6 according to the Bayesian theorem, and obtaining an expected value of the calculation center parameter by the multi-center expected integral formula expressed by the likelihood function and prior probability.

$\begin{matrix} {{\overset{\hat{}}{\gamma}}_{ig} = {\frac{1}{n_{i}}{\sum_{j}Z_{ijg}}}} & (3) \end{matrix}$ $\begin{matrix} {{\overset{\hat{}}{\delta}}_{ig}^{2} = {\frac{1}{n_{i - 1}}{\sum_{j}\left( {Z_{ijg} - {\overset{\hat{}}{\gamma}}_{ig}} \right)^{2}}}} & (4) \end{matrix}$ $\begin{matrix} {{E\left\lbrack \gamma_{ig} \right\rbrack} = {\int{\gamma_{ig}{\pi\left( {Z_{ig},\gamma_{ig},\delta_{ig}^{2}} \right)}{d\left( {\gamma_{ig},\delta_{ig}^{2}} \right)}}}} & (5) \end{matrix}$ $\begin{matrix} {{E\left\lbrack \gamma_{ig} \right\rbrack} = {\frac{1}{C\left( Z_{ig} \right)}{\int{\gamma_{ig}{L\left( {Z_{ig}{❘{\gamma_{ig},\delta_{ig}^{2}}}} \right)}{\pi\left( {\gamma_{ig},\delta_{ig}^{2}} \right)}{d\left( {\gamma_{ig},\delta_{ig}^{2}} \right)}}}}} & (6) \end{matrix}$

where, n_(i) represents a sample size of a center i; Z_(ijg) represents the standardized features, in which i represents a center, j represents a case, and g represents a feature, C(Z_(ig))=∫L(Z_(ig)|γ_(ig), δ_(ig) ²)π(γ_(ig), δ_(ig) ²)d(γ_(ig), δ_(ig) ²) L(Z_(ig)|γ_(ig), δ_(ig))=Π_(j)φ(Z_(ig), γ_(ig), δ_(ig) ², φ(Z_(ig), γ_(ig), δ_(ig) ²) is a probability density function for evaluating a certain feature Z_(ijg) about a random variable N(γ_(ig), δ_(ig) ²) that conforms to the normal distribution;

(3.2) calculating the multi-center effect parameter integral of formula 6 by a Monte Carlo integration method, defining w_(ig″)=L(Z_(ig)|{circumflex over (γ)}_(ig″), {circumflex over (δ)}_(ig″)), so that an additive multi-center effect parameter γ_(ig)*, can be estimated; using a same process to estimate a multiplicative multi-center effect δ_(ig)*; a numerical integration method-Monte Carlo integration is used to estimate the integration result of formula 6, which avoids using an analytical method to calculate the integration.

$\begin{matrix} {\gamma_{ig}^{*} = {{\hat{E}\left( \gamma_{ig} \right)} = \frac{\sum_{g^{''}}{w_{{ig}^{''}}{\overset{\hat{}}{\gamma}}_{{ig}^{''}}^{2}}}{\sum_{g^{''}}w_{{ig}^{''}}}}} & (7) \end{matrix}$ $\begin{matrix} {\delta_{ig}^{*} = {{\hat{E}\left( \delta_{ig} \right)} = \frac{\sum_{g^{''}}{w_{{ig}^{''}}{\overset{\hat{}}{\delta}}_{{ig}^{''}}^{2}}}{\sum_{g^{''}}w_{{ig}^{''}}}}} & (8) \end{matrix}$

(3.3) substituting the multi-center effect parameters γ*_(Bg) and δ*_(Bg) estimated in step (3.2) into the location-scale model to obtain the feature Y*_(Bjg) of the test center B after the multi-center effect compensation;

$\begin{matrix} {{Y_{{Bjg}^{=}}^{*}\frac{{\overset{\hat{}}{\sigma}}_{g}}{{\hat{\delta}}_{Bg}^{*}}\left( {Z_{Bjg} - {\overset{\hat{}}{\gamma}}_{Bg}^{*}} \right)} + {\overset{\hat{}}{\alpha}}_{g} + {X{{\overset{\hat{}}{\beta}}_{g}.}}} & (9) \end{matrix}$

Examples

Taking the deep convolution neural network ResNet-152 as an example, how to apply this method for multi-center effect compensation to the deep convolution neural network will be explained.

The network structure of ResNet-152 is shown in Table 1:

TABLE 1 Resnet-152 Network Structure Layer name Detail Convolution layer 1 7 * 7, stride 2 Convolution layer 2 3 * 3 max pool, stride 2 $\begin{bmatrix} {{1*1},} & 64 \\ {{3*3},} & 64 \\ {{1*1},} & 256 \end{bmatrix}*3$ Convolution layer 3 $\begin{bmatrix} {{1*1},} & 128 \\ {{3*3},} & 128 \\ {{1*1},} & 512 \end{bmatrix}*8$ Convolution layer 4 $\begin{bmatrix} {{1*1},} & 256 \\ {{3*3},} & 256 \\ {{1*1},} & 1024 \end{bmatrix}*36$ Convolution layer 5 $\begin{bmatrix} {{1*1},} & 512 \\ {{3*3},} & 512 \\ {{1*1},} & 2048 \end{bmatrix}*3$ Mean value pooling layer + 1000-dimensional fully connected layer + softmax

In an actual implementation process, a training center A used the PET/CT image data of this center to train an intelligent auxiliary diagnosis model of benign and malignant tumors based on ResNet-152. Now, this intelligent auxiliary diagnosis model was used in a test center B, a grass-roots hospital. However, due to the statistical differences between the imaging protocols used by the training center A and the test center B, if the model was directly used in the test center B, the diagnostic efficiency would be reduced to some extent. Therefore, the above problems were solved by the method for multi-center effect compensation of the present disclosure. The implementation steps were as follows:

(1) All the image data used for training in the training center A were input into ResNet-152 after training, and a feature map of the last layer of a convolution layer 5 was extracted and saved as features.

(2) All available image data of THE test center B were input into the trained ResNet-152, and the feature map of the last layer of the convolution layer 5 was extracted and saved as features.

(3) Based on the features of the two groups obtained in step (1) and step (2), the additive and multiplicative multi-center effect parameters were estimated by using the method for multi-center effect compensation proposed by the present disclosure.

(4) By using the multi-center effect parameters obtained in step (3), the features of the test center B were compensated for a multi-center effect.

(5) The compensated features of the test center B were returned to the second half of ResNet-152, that is, the mean value pooling layer behind the convolution layer 5, the 1000-dimensional fully connected layer and the softmax layer, so as to obtain the classification results of benign and malignant tumors.

After the estimated parameters of multi-center effect were saved, the steps (1) and (3) could be omitted in the subsequent diagnosis application as follow:

(1) The image data to be tested in the test center B were input into the diagnosis model ResNet-152, and the feature map of the last layer of the convolution layer 5 was extracted as features.

(2) By using the obtained additive and multiplicative multi-center effect parameters, the features obtained in step (1) were compensated for a multi-center effect.

(3) The compensated features of test center B were returned to the second half of ResNet-152, that is, the mean value pooling layer behind convolution layer 5, the 1000-dimensional fully connected layer and the softmax layer, so as to obtain the classification results of benign and malignant tumors.

Due to the multi-center effect compensation, the data of the test center B is applied to the diagnosis model of the training center A, which has a similar diagnosis efficiency to that of the training center A, and indirectly improves the generalization ability of the model.

It should be noted that when the data compression apparatus provided in the foregoing embodiment performs data compression, division into the foregoing functional modules is used only as an example for description. In an actual application, the foregoing functions can be allocated to and implemented by different functional modules based on a requirement, that is, an inner structure of the apparatus is divided into different functional modules, to implement all or some of the functions described above. For details about a specific implementation process, refer to the method embodiment. Details are not described herein again.

All or some of the foregoing embodiments may be implemented by using software, hardware, firmware, or any combination thereof. When the software is used for implementation, all or some of the embodiments may be implemented in a form of a computer program product. The computer program product includes one or more computer instructions. When the computer program instructions are loaded and executed on a server or a terminal, all or some of the procedures or functions according to the embodiments of this application are generated. The computer instructions may be stored in a computer-readable storage medium or may be transmitted from a computer-readable storage medium to another computer-readable storage medium. For example, the computer instructions may be transmitted from a web site, computer, server, or data center to another web site, computer, server, or data center in a wired (for example, a coaxial optical cable, an optical fiber, or a digital subscriber line) or wireless (for example, infrared, radio, or microwave) manner. The computer-readable storage medium may be any usable medium accessible by a server or a terminal, or a data storage device, such as a server or a data center, integrating one or more usable media. The usable medium may be a magnetic medium (for example, a floppy disk, a hard disk, or a magnetic tape), an optical medium (for example, a digital video disk (DVD)), or a semiconductor medium (for example, a solid-state drive).

The steps of the method or algorithm described combined with the embodiments of the present disclosure may be implemented in a hardware manner, or may be implemented in a manner in which a processor executes software instructions. The software instructions may consist of corresponding software modules, and the software modules can be stored in Random Access Memory (RAM), flash memory, Read Only Memory (ROM), Erasable Programmable ROM (EPROM), Electrically Erasable Programmable Read-Only Memory (EEPROM), registers, hard disks, removable hard disks, CD-ROMs or any other forms of storage media well-known in the art. An exemplary storage medium is coupled to the processor, such that the processor can read information from, and write information to, the storage medium. The storage medium can also be an integral part of the processor. The processor and storage medium may reside in an Application Specific Integrated Circuit (ASIC). Alternatively, the ASIC may be located in a node device, such as the processing node described above. In addition, the processor and storage medium may also exist in the node device as discrete components. 

What is claimed is:
 1. A method for multi-center effect compensation based on a PET/CT intelligent diagnosis system, comprising the following steps: importing PET/CT data of a training center A into a deep convolution neural network model of the training center A, obtaining a feature map of all the PET/CT data in the training center A, calculating an average value {circumflex over (α)}_(g), a covariant regression coefficient {circumflex over (β)}_(g) and a variance {circumflex over (σ)}_(g) of all features Y_(Ajg) in the feature map of all the PET/CT data in the training center A, and carrying out Z-Score standardization thereon to obtain a standardized feature Z_(Ajg): $\begin{matrix} {Z_{Ajg} = \frac{Y_{Ajg} - {\overset{\hat{}}{\alpha}}_{9} - {X{\overset{\hat{}}{\beta}}_{g}}}{{\overset{\hat{}}{\sigma}}_{g}}} & (1) \end{matrix}$ where Z_(Ajg) represents a feature g extracted from a case j in the training center A, and X is a design matrix of a relevant covariates; importing PET/CT data in a test center B into the deep convolution neural network model of the training center A, obtaining a feature map of all the PET/CT data in the test center B, and carrying out Z-Score standardization on the features Y_(Bjg) in the feature map by using the average value and the variance of the training center A, and obtaining a standardized feature Z_(Bjg): $\begin{matrix} {{Z_{Bjg} = \frac{Y_{Bjg} - {\overset{\hat{}}{\alpha}}_{g} - {X{\overset{\hat{}}{\beta}}_{g}}}{{\overset{\hat{}}{\sigma}}_{g}}};} & (2) \end{matrix}$ estimating multi-center effect parameters γ*_(ig), and δ*_(ig), according to the standardized features of the training center A and the test center B, and carrying out multi-center effect compensation for the feature of the test center B, so as to obtain a feature Y*_(ijg) of the test center B after the multi-center effect compensation, which comprise following sub-steps: the standardized feature Z_(Ajg) and the standardized feature Z_(Bjg) conforming to a normal distribution, and calculating an average value {circumflex over (γ)}_(ig) and a standard deviation {circumflex over (δ)}_(ig) ² of the standardized features; obtaining an expected value E[γ_(ig)] about a multi-center effect parameter γ_(ig) by integration from a posterior distribution π(Z_(ig), γ_(ig), δ_(ig) ²) of standardized data Z_(ig); transforming the integral by a likelihood function according to the Bayesian theorem; $\begin{matrix} {{\overset{\hat{}}{\gamma}}_{ig} = {\frac{1}{n_{i}}{\sum_{j}Z_{ijg}}}} & (3) \end{matrix}$ $\begin{matrix} {{\overset{\hat{}}{\delta}}_{ig}^{2} = {\frac{1}{n_{i - 1}}{\sum_{j}\left( {Z_{ijg} - {\overset{\hat{}}{\gamma}}_{ig}} \right)^{2}}}} & (4) \end{matrix}$ $\begin{matrix} {{E\left\lbrack \gamma_{ig} \right\rbrack} = {\int{\gamma_{ig}{\pi\left( {Z_{ig},\gamma_{ig},\delta_{ig}^{2}} \right)}{d\left( {\gamma_{ig},\delta_{ig}^{2}} \right)}}}} & (5) \end{matrix}$ $\begin{matrix} {{E\left\lbrack \gamma_{ig} \right\rbrack} = {\frac{1}{C\left( Z_{ig} \right)}{\int{\gamma_{ig}{L\left( {Z_{ig}{❘{\gamma_{ig},\delta_{ig}^{2}}}} \right)}{\pi\left( {\gamma_{ig},\delta_{ig}^{2}} \right)}{d\left( {\gamma_{ig},\delta_{ig}^{2}} \right)}}}}} & (6) \end{matrix}$ where n_(i) represents a sample size of a center i; Z_(ijg) represents the standardized features, in which i represents a center, j represents a case, and g represents a feature, C(Z_(ig))=∫L(Z_(ig)|γ_(ig), δ_(ig) ²)π(γ_(ig), δ_(ig) ²)d(γ_(ig), δ_(ig) ²), L(Z_(ig)|γ_(ig), δ_(ig) ²)=Π_(j)φ(Z_(ijg), γ_(ig), δ_(ig) ²), φ(Z_(ijg), γ_(ig), δ_(ig) ²) is a probability density function for evaluating a certain feature Z_(ijg) about a random variable N(γ_(ig), δ_(ig) ²) that conforms to the normal distribution; calculating an integral about the multi-center effect parameters by a Monte Carlo integral method, defining w_(ig″)=L(Z_(ig)|{circumflex over (γ)}_(ig″), {circumflex over (δ)}_(ig″)), so that an additive multi-center effect parameter γ*_(ig) can be estimated; using a same process to estimate a multiplicative multi-center effect δ*_(ig); $\begin{matrix} {\gamma_{ig}^{*} = {{\hat{E}\left( \gamma_{ig} \right)} = \frac{\sum_{g^{''}}{w_{{ig}^{''}}{\overset{\hat{}}{\gamma}}_{{ig}^{''}}^{2}}}{\sum_{g^{''}}w_{{ig}^{''}}}}} & (7) \end{matrix}$ $\begin{matrix} {\delta_{ig}^{*} = {{\hat{E}\left( \delta_{ig} \right)} = \frac{\sum_{g^{''}}{w_{{ig}^{''}}{\overset{\hat{}}{\delta}}_{{ig}^{''}}^{2}}}{\sum_{g^{''}}w_{{ig}^{''}}}}} & (8) \end{matrix}$ substituting the multi-center effect parameters γ*_(Bg) and δ*_(Bg) into the location-scale model to obtain the feature Y*_(Bjg) of the test center B after the multi-center effect compensation; $\begin{matrix} {{Y_{{Bjg}^{=}}^{*}\frac{{\overset{\hat{}}{\sigma}}_{g}}{{\hat{\delta}}_{Bg}^{*}}\left( {Z_{Bjg} - {\overset{\hat{}}{\gamma}}_{Bg}^{*}} \right)} + {\overset{\hat{}}{\alpha}}_{g} + {X{{\overset{\hat{}}{\beta}}_{g}.}}} & (9) \end{matrix}$ 